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Abstract ofFR2775094 

The operation of delumping is based on the fact that the oil reservoir has only two non aqueous phases 
(liquid and gas) , allowing fluid regrouping and detailed decomposition calculations to be carried out. 
The oil reservoir is split into a grid system and aggregate fluid compositional simulation using 
components and pseudo components carried out. For each grid section, modeling is carried out on the 
vaporized fraction with the regrouped fluid under constant equilibrium. For each grid step, a molar 
quantity is produced and evaluation of the detailed fluid composition carried out. 
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METHODE DE SIMULATION POUR PREDIRE EN FONCTION DU TEMPS UNE COMPOSITION DETAILLEE 
D'UN FLUIDE PRODUIT PAR UN RESERVOIR. 



Methode de simulation pour predire en fonction du 
temps une composition detaillee a Q composants et/ ou 
pseudo-composants d'un fluide produit par un reservoir. 
Elle est du type consistant a: 

- representer le reservoir sous la forme d'un reseau de 
mailles G) dont chacune constitue un volume elementaire 
rempli de fluide, 

- realiser, de fagon connue en soi, une simulation com- 
positionnelle du fluide regroupe a N composants et pseudo- 
composants (i), N etant inferieur a Q, et elle est caracterisee 
en ce qu'elle consiste en outre a: 

e) determiner, a chaque pas le temps (m) et pour cha- 
que maille (j) la composition des phases liquide et vapeur du 
fluide detaille de la composition Q du fluide, a partir des va- 
leurs de la fraction vaporisee (Oj" 1 ) et des constantes d'equi- 
libre (kij m ) du fluide regroupe, 

f) evaluer pour chaque maille, au pas de temps (m+1 ), la 
quantite molaire de chacun des Q composants et/ ou pseu- 
do-composants du fluide detaille a partir des valeurs corres- 
pondantes au pas de temps (m), des debits des phases du 
fluide regroupe et des compositions des phases liquide et 
vapeur du fluide detaille determinees a retape e), et a 

g) evaluer, pour chaque puits de production, la compo- 



sition detaillee du fluide produit entre des instants t et t' cor- 
respondant aux pas de temps m et m+1 , a partir des debits 
des phases du fluide regroupe et des compositions des pha- 
ses liquide et vapeur du fluide detaille determinees a I'etape 
e). 
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Methode de simulation pour predire en fonction du temps 
une composition detaillee d'un fluide produit par un reservoir 



5 

La presente invention conceme une methode de simulation pour 
predire, en fonction du temps, la composition detaillee d'un fluide produit 
par un reservoir et plus particulierement la composition detaiUee d'un fluide 
contenu dans et produit par un gisement petrolier dans lequel sont implantes 

10 unou plusieurs puits de production. 

La simulation compositionnelle d'un gisement petrolier est 
couramment utilisee pour fournir les profils previsionnels de production du 
gisement qui permettent notamment de determiner le schema de production 
le mieux adapte a ce gisement. 

15 La simulation compositionnelle d'un gisement petrolier est mise 

en oeuvre non pas en utilisant une description reelle du fluide mais en 
utilisant un fluide modelise par un nombre de composants plus reduit que le 
nombre de composants du fluide reel. En effet, le nombre de composants du 
fluide de gisement etant relativement grand, la moderation avec tous les 

20 composants conduirait a des temps de calcul trop importants. Ce temps de 
calcul prohibitif a amend les specialistes a regrouper les composants purs en 
pseudo-composants, par exemple des pseudo-composants regroupant Tazote 
et le methane, des pseudo-composants regroupant les hydrocarbures en C 3 
et C 4 , . . . etc. . et a realiser la simulation de gisement sur un nombre reduit 

25 N de composants purs et de pseudo-composants. Une telle composition 
reduite est designee par l f expression "composition regroupee du fluide" ou 
"fluide regroupe". Generalement N represente 5 a 8 composants et pseudo- 
composants, ce qui est considere comme etant suffisant pour bien 
representer le comportement du fluide de gisement aux conditions de fond. 

30 Avant d'effectuer la simulation compositionnelle de gisement, on 

represente le reservoir ou gisement sous la forme d'un reseau de mailles 
dont chacune constitue un volume elementaire dudit reservoir. Le nombre 
de mailles peut atteindre plusieurs milliers et chaque maille presente des 
propnetes qui lui sont propres, comme la geometrie, la porosite, la 
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permeabilite, etc.. . De plus, a certaines mailles peut correspondre au moins 
un puits d'injection ou de production qui est implante dans ledit reservoir. 

La simulation de gisement permet de calculer pour chaque maille 
un certain nombre de variables principales du fluide de gisement, lesdites 
5 variables pouvant etre la quantite (nombre de moles) de fluide, la compo- 
sition ou fraction molaire de chaque composant et pseudo-composant du 
fluide regroupe et la pression regnant dans chaque maille. Ces variables 
sont connues a l'instant t=0 (debut de Texploitation du gisement) et sont 
ensuite calculees par la simulation de gisement a chaque instant t ou pas de 

10 temps m. Pour chaque maille, il est possible, a partir de ces variables 
principales, de calculer toute autre propriete du fluide present dans la 
maille, comme le nombre de phases, la composition de chaque phase, etc. 
Pour chaque maille a laquelle est associe un puits de production, il est 
possible, a partir de ces variables principales et des contraintes de 

15 production imposees, de calculer egalement le debit de production et la 
composition du fluide produit par ledit puits. 

Les proprietes thermodynamiques du fluide peuvent etre 
calculees en utilisant une des equations d'etat bien connues des specialistes 
et qui ne seront done pas rappelees. 

20 La modelisation a N composants, si elle est suffisante pour 

representer le comportement du fluide de gisement aux conditions de fond, 
n'est plus appropriee pour simuler le comportement du fluide dans les 
installations pour son exploitation en surface, laquelle exploitation necessite 
de connaitre une composition plus detaillee a Q composants et/ou pseudo- 

25 composants du fluide produit par le gisement, Q etant superieur a N et par 
exemple de l'ordre de 20 a 30. 

Jusqu'ici, on obtenait cette composition detaillee a Q composants 
et/ou pseudo-composants a partir de la composition regroupee N en 
supposant que la composition des pseudo-composants restait constante dans 

30 le temps. 

Une telle maniere de proceder est source d'erreurs car la 
composition des pseudo-composants varie dans le temps, notamment en 
fonction de la pression du gisement ou suite a rinjection d'un gaz dans le 
gisement. 
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De ce fait, la composition detaillee obtenue reste tres approxima- 
tive et ne permet pas de predire de maniere satisfaisante le comportement du 
fluide dans les conditions Sexploitation en surface bien des annees apres le 
debut de la production. 
5 La presente invention a pour but de proposer une methode de 

simulation pour predire la composition detaillee du fluide produit par un 
reservoir en tenant compte des parametres du reservoir et qui soit beaucoup 
plus precise. 

L'invention a pour objet une methode du type consistant a : 
10 a) representer le reservoir sous la forme d f un reseau de mailles (j) 

dont chacune constitue un volume elementaire rempli de fluide, 

b) defmir le fluide par une modelisation regroupee a N 
composants et pseudo-composants (i) et a determiner une equation d'etat 
decrivant le fluide dans cette modelisation regroupee, 
15 c) definir egalement le fluide par une modelisation detaillee a Q 

composants et/ou pseudo-composants, Q etant superieur a N, et a deter- 
miner une autre equation d'etat decrivant le fluide dans cette modelisation 
detaillee, 

d) realiser, de fa9on connue en soi, une simulation composi- 
20 tionnelle du fluide regroupe a N composants et pseudo-composants (i), 

ladite simulation compositionnelle permettant de calculer au moins pour 
chaque maille (j) et a des pas de temps consecutifs (m, m+1, ...) la fraction 
vaporisee (0™), les constantes d'equilibre liquide-vapeur (k™ ) de chaque 

composant (i), les debits d'injection ou de production (s™) et pour chaque 
25 paire de mailles (j, h) les debits des phases liquide (u^ h ) et vapeur 

(^h) du 

fluide a N composants et pseudo-composants, et elle est caracterisee en ce 
qu'elle consiste en outre a: 

e) determiner, a chaque pas de temps (m) et pour chaque maille 
(j) la composition des phases liquide et vapeur du fluide detaille de la 

30 composition Q du fluide, a partir des valeurs de la fraction vaporisee (9™) et 

des constantes d'equilibre (kJJ 1 ) du fluide regroupe, 

f) evaluer pour chaque maille, au pas de temps (m+1), la quantite 
molaire de chacun des Q composants et/ou pseudo-composants du fluide 
detaille a partir des valeurs correspondantes au pas de temps (m), des debits 
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des phases du fluide regroupe et des compositions des phases liquide et 
vapeur du fluide detaille determinees a l'etape e), et a 

g) evaluer, pour chaque puits de production, la composition 
detaillee du fluide produit entre des instants t et t correspondant aux pas de 
5 temps m et m+1, a partir des debits des phases du fluide regroupe et des 
compositions des phases liquide et vapeur du fluide detaille determinees a 
l'etape e). 

Selon une autre caracteristique de la presente invention, les 
etapes e, f et g sont mises en oeuvre en meme temps que la realisation de la 

10 simulation compositionneile du fluide regroupe. 

Selon une autre caracteristique de la presente invention, les 
resultats de l'etape d sont memorises puis utilises ulterieurement pour la 
mise en oeuvre des etapes e, f et g. 

La methode selon la presente invention est mise en oeuvre pour 

15 un reservoir constitue, par exemple, par un gisement petrolier. De maniere 
connue, on realise une simulation de gisement. A cet effet, on represente le 
gisement sous la forme d'un reseau maille, certaines des mailles ou un 
groupe de mailles etant associes a un des puits de production implantes dans 
le gisement petrolier a exploiter. 

20 Pour reduire le temps de calcul, la simulation compositionneile 

du gisement est effectuee sur un nombre limite N de composants et pseudo- 
composants, par exemple de 5 a 8, definis de la maniere indiquee 
precedemment, ces composants et pseudo-composants etant selectionnes en 
fonction de la nature du gisement. 

25 Le fluide regroupe a N composants et pseudo-composants est 

decrit par une equation d'etat qui pourrait etre par exemple celle de PENG- 
ROBINSON adaptee a ladite composition regroupee. 

Au debut de Texploitation du gisement, on connait, par des 
mesures effectuees au prealable, pour chaque maille j du reseau, les 

30 variables principals a l'instant t=0 ou pas de temps m=0. Ces variables 
principales sont : 
la quantite de fluide : f ™ 

la fraction molaire du composant i : z™ avec 1 < i < N 
la pression : p" 
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Pour passer du pas de temps m au pas de temps m+1, la 
simulation compositionnelle comprend plusieurs etapes detaillees ci-apres. 

Dans une premiere etape et pour le pas de temps m, un calcul 
d'equilibre liquide-vapeur (flash) est effectue sur le fluide regroupe, pour 
5 chaque maille j afin de determiner : 

la fraction vaporisee 

la fraction molaire de chaque composant i dans la phase huile (x™ avec i = 1...n) et 
la fraction molaire dudit composant i dans la phase gaz (y™ avec i = 1...n). 

Ces diverses fractions sont determinees avec le systeme 
10 d'equations suivant : 

4>i(xJ,....x- ) et * i (yJ,...-yS i ) avec 1 < i < N 

(i-e:)xj+e- y ?=2? 

et 

N N 

1=1 1=1 

15 

dans lesquelles : 

<t>{ est la fugacite du composant i. 

Selon une premiere caracteristique de I'invention, on memorise la 

y m 

fraction vaporisee 0* et les constantes d'equilibre k™=— ^ avec i variant 

20 entre 1 et N, afin de les utiliser dans une autre etape denommee 
"delumping" qui sera explicitee ulterieurement 

Dans une deuxieme etape et pour chaque maille, on evalue 
certaines des proprietes de chaque phase liquide et gaz, a savoir masse 
volumique (p™ et p™), viscosite, saturation et permeabilite relative. 

25 Dans une troisieme etape, pour chaque maille j et au pas de 

temps m, on calcule les coefficients de l'equation pour la pression. 

Les equations pour la pression d'une maille j a Tinstant t+At ou au 
pas de temps m+1, sont : 
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Pi "Pi )- 



At 



Pj 



m+1 



Pour resoudre ces equations, les coefficients a calculer sont : 
5 - la compressibilite totale C"j (fluide + roche), 

- les transmissivites generalisees 

T5 = ZaS , (X 0 x i +X,y i )" t jh 

1=1 

N 

10 - le tenne source : s™ = £ s ™ 

i=i 

Dans ces equations : 
h sont les mailles voisines de la maille j, h e J(j), 

a™ est le volume molaire partiel du composant i dans la maille j au pas de 
temps m, et A™ sont les mobilites des phases huile (o) et gaz (g), 
15 t jh est la transmissivite entre les centres des mailles j et k, 1 etant la maille 

amont. 

Dans une quatrieme etape, on resout pour toutes les mailles les 
equations pour la pression a Finstant t+At c'est-a-dire au pas de temps m+1. 
Ces equations forment un systeme lineaire du type : 

20 

A. p = B 

a a; ~ 

ou p = (p^ 1 -p^* 1 )\ J etant le nombre de mailles du reseau. 

Ces equations sont resolues par un procede iteratif standard. 
25 Dans une cinquieme etape, on calcule au pas de temps m+1, les 

debits des phases entre chaque paire de mailles adjacentes et pour chaque 
puits associe a partir, par exemple, de la loi de DARC Y. 

^^-Po-^-tjpr-pr 1 ) 
ui=-p^(pr , -pr) 

30 
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1 etant la maille amont avec l=j ou h. 

Selon une autre caracteristique de 1'invention, on stocke ou 
memorise les debits des phases qui seront utilises dans le delumping. 

Dans une sixieme etape, on evalue la quantite f et la composition 
z du fluide regroupe, dans chaque maille j et au pas de temps m+ 1 a partir 
des valeurs au pas de temps m, des debits des phases evalues a la cinquieme 
etape et des compositions des phases calculees a Tetape 1. Pour ce faire, on 
applique les lois de conservation defrnie par les equations suivantes : 

f-'z^ 1 = fx - At S( x u^h + y>S,)-sj 

Par un processus iteratif, on repete les etapes ci-dessus pour 
calculer les variables principales aux temps t+2At, t+3At, etc... ou encore 
aux pas de temps m+1, m+2, m+3, etc., l'increment At pouvant etre 
constant ou variable. 

Selon la presente invention, on defmit une modelisation detaillee 
a Q composants et/ou pseudo-composants ainsi qu'une autre equation d'etat 
decrivant le fluide ainsi modelise. 

Pour cela, on choisit prealablement les composants et/ou pseudo- 
composants dont on veut connaitre les caracteristiques et qui sont 
necessaires pour representer le fluide de production dans les conditions de 
surface. Le nombre Q qui est superieur au nombre N de la composition 
regroupee est generalement de Tordre de 20 a 30. On connait la composition 
du fluide detaille a Q composants et/ou pseudo-composants au pas de temps 
m=0. 

Les donnees stockees aux premiere et cinquieme etapes sont 
utilisees pour les etendre aux Q composants et/ou pseudo-composants de la 
composition detaillee. 

Comme Tecoxilement peut s'effectuer dans n'importe quelle 
direction entre les mailles voisines, chaque maille ne peut pas etre traitee 
independamment. C'est pourquoi, a chaque pas de temps, toutes les mailles 
sont traitees ensemble. 
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Pour realiser le delumping, on fait ttiypothese qu'a chaque pas de 

temps m, 

- la fraction molaire de la phase vaporisee dans chaque maille est 
independante du nombre de constituants du fluide, 

- le nombre de moles de liquide et de vapeur s'ecoulant entre les mailles et 
dans ou hors des puits est independant du nombre de constituants du 
fluide, 

U3. = ujrt us, = u™ 
S5 = s5 et S£ = s£ 

Les lettres en majuscule sont utilisees pour representer le fluide 
detaille et les lettres en minuscule sont utilisees pour representer le fluide 
regroupe. 

Le calcul des constantes d'equilibre peut etre effectue a partir de 
15 1'expression : 

Ln[K,] = AC 0 + £ AC^ ou AC n = C L Q - C n v 

comme indique par C. LEIBOVICI dans son article "A Consistent 
Procedure for Pseudo-Component Delumping" paru dans FLUID PHASE 
20 EQUILIBRIA, 117 (1996), 225-232. 

Pour une equation d'etat a deux parametres, comme celle de 
PENG ROBINSON, l'equation generalisee de LEIBOVICI peut etre ecrite 
pour le fluide detaille de la maniere suivante : 

25 Ln(K I 7) = T 1 --fP^-fyrb I 

ou aj et b\ sont les parametres de Tequation d'etat pour le composant I et ou 
les parametres r^&J ety™ sont calcules a partir des constantes d'equilibre 

pour le fluide regroupe (k™) en minimisant la fonction : 



5 
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De ce fait, lorsqu'on connait la composition detaillee z™ dans la 
maille j au pas de temps in, on peut estimer les compositions detaillees des 
phases liquide et vapeur dans la meme maille j et pour le meme pas de 
temps m. 

5 

am K m Z m 

X?= £ r Y"= M r avecl<I<Z 

I etant le composant et/ou pseudo-composant de la composition 
detaillee. Ces compositions peuvent etre normalisees afm que leur somme 
10 soit egale a 1 si necessaire. 

Lorsqu'on connait les compositions des phases detaillees pour 
toutes les mailles du reseau au pas de temps m, il est possible d'estimer les 
compositions detaillees au pas de temps consecutif m+1 a partir des 

equations ci-apres : 



15 



z?ft-(y5S3+x5s;)- Z (yTO + xTO) 



zr'= avecI=l,...,Q 



ou 



20 F? +i =Fr-(s;+s;)- zte+u^h) 

1 heJ(j) 

dans lesquelles j -j pour un debit positif de la maille j vers la maille h ou 
dans le puits et j-h pour un debit negatif, j' correspondant au fluide injecte 
pour les puits ^injection ou S sera negatif. 

Grace a la presente invention, on obtient une composition 

25 detaillee a Q composants et/ou pseudo-composants qui permet de mieux 
definir les profils de production previsionnels et ce, avec une precision 
accrue comparee aux methodes anterieures. Cette plus grande precision des 
profils de production previsionnels permet un choix de schema de 
developpement (type et taille d'installations de surface, nombre de puits, 

30 etc..) plus fiable, avec pour consequence une augmentation de la rentabilite 
et une diminution du risque economique. 
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Dans ce qui precede, certaines valeurs obtenues au cours de la 
simulation de gisement sont stockees ou memorisees pour etre utilisees 
ulterieurement dans le delumping. II va de soi que ces memes valeurs 
pourraient etre utilisees directement si on decide d'effectuer le delumping en 
meme temps que ladite simulation de gisement. 
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REVINDICATIONS 



5 1. Methode de simulation pour predire, en fonction du temps, la 

composition detaillee d'un fluide produit par un reservoir, du type consistant 

a : 

a) representer le reservoir sous la forme d'un reseau de mailles (j) 
dont chacune constitue un volume elementaire rempli de fluide, 
10 b) definir le fluide par une modelisation regroupee a N compo- 

sants et pseudo-composants (i) et a determiner une equation d'etat decrivant 
le fluide dans cette modelisation regroupee, 

c) definir egalement le fluide par une modelisation detaillee a Q 
composants et/ou pseudo-composants, Q etant superieur a N, et a deter- 

15 miner une autre equation d'etat decrivant le fluide dans cette modelisation 
detaillee, 

d) realiser, de fa$on connue en soi, une simulation composi- 
tionnelle du fluide regroupe a N composants et pseudo-composants (i), 
ladite simulation compositionnelle permettant de calculer au moins pour 

20 chaque maille (j) et a des pas de temps consecutifs (m, m+1, ...) la fraction 
vaporise (0™), les constantes d'equilibre liquide-vapeur (k?) de chaque 
composant (i), les debits d'injection ou de production (s™) et pour chaque 
paire de mailles (j, h) les debits des phases liquide (u" h ) et vapeur 
fluide a N composants et pseudo-composants, et elle est caracterisee en ce 

25 qu'elle consiste en outre a : 

e) determiner, a chaque pas le temps (m) et pour chaque maille 
(j) la composition des phases liquide et vapeur du fluide detaille de la 
composition Q du fluide, a partir des valeurs de la fraction vaporisee (g™) et 

des constantes d'equilibre (k^) du fluide regroupe, 

30 f) evaluer pour chaque maille, au pas de temps (m+1), la quantite 

molaire de chacun des Q composants et/ou pseudo-composants du fluide 
detaille a partir des valeurs correspondantes au pas de temps (m), des debits 
des phases du fluide regroupe et des compositions des phases liquide et 
vapeur du fluide detaille determiiiees a l'etape e), et a 
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g) evaluer, pour chaque puits de production, la composition 
detaiilee du fluide produit entire des instants t et t 1 correspondant aux pas de 
temps m et m+1, a paitir des debits des phases du fluide regroupe et des 
compositions des phases liquide et vapeur du fluide detaille determinees a 
5 Tetape e). 

2. Methode selon la revendication 1, caracterisee en ce que les 
etapes e, f et g sont mises en oeuvre en meme temps que la realisation de la 
simulation compositionnelle du fluide regroupe. 

3. Methode selon la revendication 1, caracterisee en ce que les 
10 resultats de Tetape d sont memorises puis utilises ulterieurement pour la 

mise en oeuvre des etapes e, f et g. 
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